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Ch , 

' Two-level boson systems displaying a quantum phase transition from a spherical (symmetric) to 

a deformed (broken) phase are studied. A formalism to diagonalize Hamiltonians with 0{2L + 1) 
symmetry for large number of bosons is worked out. Analytical results beyond the simple mean-field 
treatment are deduced by using the continuous unitary transformations technique. In this scheme, a 
1 /N expansion for different observables is proposed and allows one to compute the finite-size scaling 
. exponents at the critical point. Analytical and numerical results are compared and reveal the power 

O ' of the present approach to compute the finite-size corrections in such a context. 
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I. INTRODUCTION 



The study of two-level systems has been a topic of interest since the first steps in the development of quantum 
''O ■ mechanics. The main advantage of these models is that they can be numerically diagonalized for very large dimensions 
Ch ' and, at the same time, they can model realistic quantum many-body systems. Typical examples are the Jaynes- 
O \ Cummings model of quantum optics the Vibron Model (VM) of quantum chemistry the two-level pairing 
model in condensed matter Q and in nuclear physics 0|, the Lipkin-Meshkov-Glick model (LMG) 0,0,0 and the 
\ Interacting Boson Model (IBM) jsj of nuclear structure. While some of these models describe two-level fermion 
. systems, the model Hamiltonian can always be written in terms of SU{2) pseudo-spin operators. Subsequently, the 
spin Hamiltonian can be expressed in terms of bosons using either the finite Schwinger representation or the infinite 



Holstein-Primakoff representation of the SU{2) algebra. An example is the LMG model which has recently been newly 



revived as a model of quantum spins with l ong -range interactions to investigate the relationship between entanglement 
and quantum phase transitions (QPTs) P. llCllTlllo.ll3LIT^IlH . In its boson representation, it has also been recently 
' used as a simplified model to describe the Josephson effect between two Bose-Einstein condensates [T^ . 
lO In this work, we focus on finite two-level boson Hamiltonians, having the common feature that the lower level 

• is always a scalar L = boson, hereafter written as s boson. The upper level can have different multipolarities, 
generically noted as L whose value defines a particular model. The LMG model in the Schwinger representation has a 
second scalar {L = 0) boson for the upper level. A dipolar (L = 1) boson leads to the VM and a quadrupolar {L = 2) 
boson corresponds to the IBM. Higher angular momentum bosons can lead new models, like for example a model of 
octupole vibrations in terms of s and octupolar {L — 3) bosons. A schematic representation of the model is shown in 
g Fig. □ 
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FIG. 1: Schematic representation of the two-level boson model studied in this paper. 
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All these two-level boson models are governed by an algebraic structure which is constructed out of all bilinear 
combinations of creation and annihilation boson operators which generate the algebra of U{2L + 2). One of the 
main features of these models is that one can construct dynamical symmetries in which the Hamiltonian can be 
written in terms of Casimir invariants of a nested chain of subalgebras of U{2L + 2). In these particular cases, the 
problem is analytically solvable providing a powerful tool to check approximate methods and a reference for more 
detailed calculations. In addition, if the Hamiltonian is written as a linear combination of Casimir invariants of the 
subalgebras U{2L + 1) and 0{2L + 2) the model is still quantum integrable but requires then to solve Bethe-like 
equations numerically ,17j . The exact solution, given by Richardson almost forty years ago [18,] . is reduced to a set 
of M nonlinear coupled equations, where M is the number of boson pairs. For two-level boson models it turns out 
that the numerical diagonalization of the Hamiltonian presented below is more efficient than solving the Richardson 
equations. 

The aim of this work is to study the QPT that occurs in the two-level boson system as it evolves from the spherical 
vibrational U{2L + 1) symmetry to the deformed 0(2L-|-2) symmetry, as a function of a control parameter. Although, 
strictly speaking QPTs are defined for macroscopic systems, there is a renewed interest in studying structural changes 
in finite-size systems as the precursors of a QPT in the thermodynamic limit. Traces of these QPTs are re adily 
observed in finite systems and their properties are then correlated with the idealized thermodynamic system 
The understanding of the modifications on the characteristics of the QPT induced by finite-size effects is of crucial 
importance to extend the concept of phase transitions to finite systems. Several techniques have long been used 
to extrapolate numerical results obtained by large-scale diagonalizations or Monte Carlo calculations to the infinite 
system. Here, we focus on a somewhat inverse problem which is the finite-size corrections to the observables in two- 
level boson models like the ground-state energy, the gap, the occupation number and some electromagnetic transition 
rates. While the zeroth order in the boson number N is given by the Hartree mean-field approach for the ground 
state and the Random Phase Approximation for the excited states, going beyond this order implies the use of more 
sophisticated techniques. We make use here of the Continuous Unitary Transformations (CUTs) and give the first 
corrections to the observables in the whole U{2L + 1) to 0{2L + 2) transition. 

The structure of the paper is the following. In Sec. ^wc introduce the two-level boson models and the formalism 
for the numerical diagonalization for very large number of bosons. Section IIIII describes the mean-field treatment 
of the two-level boson models. Section Hvl is devoted to the study of the symmetric phase using CUTs. Analytical 
expressions for different orders in the 1/iV expansion of the ground-state energy, the gap, the expectation value of 
the number of L bosons in the ground state and the transition matrix element between the ground state and the first 
excited state are obtained. In Sec. 0the broken phase is analyzed, and in Sec. ^Jthe study of the critical point 
is presented from the spherical phase. In this section, we obtain the finite-size scaling exponents for the quantities 
cited above by analyzing the divergence of their expansion. In Sec. IVIII a comparison of the numerical results 
obtained using the formalism presented in Sec. |n]with the analytical CUTs results is presented. Section rVIIII is for 
summary and conclusions. Technical details concerning fiow equations can be found in appendices. 



In this Section, we present a simple algorithm for diagonalizing boson Hamiltonians with 0{2L + 1) symmetry for 
large boson numbers. The formalism is based on previous studies (2CL .21,) and is a generalization of the one presented 
recently for treating the IBM [2^ . 

We consider the following boson pairing Hamiltonian 



II. TWO-LEVEL BOSON MODELS 




(1) 



with 




(2) 



fi— — L 





(3) 



+L 




(4) 



Ljj^ creates a boson in the excited L level with projection /i, while destroys it. We have introduced above the pair 
operators with p = s or L, which will be used later on. The total number of bosons N = ns + ul is a. conserved 
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quantity. For x = 0, H can be cast into a linear combination of the quadratic Casimir operators of 0{2L + 2) and 
the corresponding subalgebras, whereas for x = 1, H is the linear Casimir operator of the U{2L + 1) algebra. Here, 
X plays the role of a control parameter, driving the system from the 0{2L + 2) deformed phase to the U{2L + 1) 
spherical phase. 

The boson pairing Hamiltonian can be studied in an elegant way by means of the noncompact SU{1, 1) algebra 
of boson pair operators. For the subspace of p bosons, where p stands generically either for s or i bosons, the SU{1, 1) 
generators are the raising operator the lowering operator K~ = {K^^ and the Cartan operator defined as 
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(5) 
(6) 
(7) 



with Dp = 2p+ 1. The three operators {A'+, K , Kp} satisfy the SU{1, 1) commutator algebra 



K^^Kp, 



= ±5p.p,Kf, 
= -26pyK°p. 



(8) 
(9) 



The Casimir operator is 



2 (^^"^P- 



D„ (D 



(10) 



The complete set of eigenstates of the pairing Hamiltonian can be constructed as a direct product of subspaces 
associated to s and L bosons. Each subspace can be written in terms of the rising operators acting on the 
corresponding subspace of unpaired p bosons 



-.K 



'^P.Vo 



0, Vp) , 



(11) 



where Vs = vl=o = 0, 1 and vl^q = 0, 1, 2, . . . The quantity Vp is known as the boson seniority for p bosons and gives 
the number of bosons of type p not coupled in pairs to zero. Note that from now on the label n means number of 
boson pairs coupled to zero angular momentum. The total number of bosons is 2fis + 2hL + Vs + VL. The normalization 
constant in (|ll|l can be obtained from the action of K~ and on the p subspacc \np = 0, i^p) 



Kp {hp 



0, 1'p 



0, 



1 

4^' 
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(12) 
(13) 



and the commutation relation 



[[K;,K;],Kj]^2R- 



(14) 



then 



0,vp 



^ \nr. 



0,i/, 



p/ ' 



(15) 



and 



(np, Vp\np, Vp) = {2hp + 2p + 2vp - 1) {fip - 1, Vp\np - 1, Vp), 



(16) 



and finally 



{np,Vp\np,Vp) 



np\{2np + 2p + 2vp-l)\\ 



2"p {2p + 2vp-l)\ 



(17) 
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Remember that the label p stands for s or L bosons and takes numerical values: for s bosons and L for L bosons. 
Once the basis for each subspace is obtained H11I17|) . the complete basis set for the pairing Hamiltonian is easily 
constructed, 



/-iris /-iriL 



{Ktr {Ktt" \f^s - 0, V,; % - 0, vl) . (18) 



We now diagonalize the Hamiltonian Q in the basis (|18(l . Note that in the following n refers to boson number 
operators but n are number of boson pairs. The relevant matrix elements for the construction of the Hamiltonian 
matrix are: 

{ns,Vs]nL,i'L\ ns \ns,i's;nL, J^l) = + Vs, (19) 
{ns,i's\nL,i'L\nL jris, i/^; fj-L, i^i) = 2fj,L + (20) 



{fis^Vs] nL,VL\KJK^ \hs, Vs;nL, vl) = fis { fis + Vs - ] , (21) 



. i^s] nL,i^L\ KJK^ \ns,Vs\nL,VL) ^ ul [ hl + L + - - ] , (22) 



{{ns - {uL + 1),vl\KIk^ \ns,Vs;nL,VL) = -\Jhs {ul + 1) (2ns + ^Vs - 1) (2nL + 2L + 2vl + 1). (23) 

It is clear that the Hamiltonian does not mix states with different boson seniority quantum numbers. Thus, the 
Hamiltonian matrix is block diagonal. In addition, within one block the matrix is tridiagonal making the diago- 
nalization simple. The different states are obtained as follows: one starts with the boson subspace containing N jl 
boson pairs coupled to zero angular momentum, Vs — vi^ — 0. The diagonalization of H in this subspace provides 
states with angular momentum zero and the first eigenstate is the ground state. Next, one goes to the block with 
one broken boson pair. This block is composed of two separate blocks since the two bosons can be one s boson and 
one L boson {vs = 1, i^l = 1) or two L bosons not coupled to zero since the coupling to zero is included in the first 
block, {vg = 0, ^'i — 2). Notice that two unpaired s bosons are not possible since they are always coupled to zero 
and consequently they are counted in the first block. For the case of the LMG model in which L = 0, the block 
j's = 0, i^L = 2 is not allowed for the same reason. The block — \-,vi^ — \ provides states with angular momentum 
L, the first of them is the first excited state of the system. The block = 0, z^l = 2 contains states with two L 
bosons not coupled to zero angular momentum, it contains angular momenta: 2L, 2L — 2, . . . , 2. Next, there is another 
block with two broken boson pairs composed of two sub-blocks: Vs — ^,i^L — S and i^s — OjVl — 4. Again, the block 
i/s = 0, i^L = 4 is absent for the LMG model. This construction continues for 3,4,..., N/2 broken boson pairs. Few 
first low-lying Hamiltonian eigenstates are depicted schematically in Fig. [21 

Direct block diagonalization of the Hamiltonian in the basis (|18|l provides observables as the ground-state energy 
or the gap and also the wave functions of the states. With the wave function of the ground state the expectation 
value of the number of L bosons in the ground state can be easily calculated. One can also calculate the transition 
probability from the ground state to the first excited state provided with the appropriate operator. The natural 
transition operator for the pairing Hamiltonian we are considering is 

TL,=Lls + {-ir.s^L^^. (24) 
The action of T^^ on the subspace Vs — O^i^L = that includes the ground state is 



r^Jn.,0;%,0)^ ^ nsiKtr {KtT' + fi^ {KtT' {K+Y' 's^lI |n, = 0, 0; n^ = 0, 0) , 



(25) 
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4L,4L-2,...,0 



2 broken boson pairs 



FIG. 2: Schematic representation of the level sequence obtained by diagonalizing the Hamiltonian for the two-level boson 
models studied in this paper. Numbers above the lines are angular momenta. 

where the state jn^ = 0,i/s = 0; = 0,1/^ = 0) = |0) is the boson vacuum. 

Then, the matrix elements of T^^ connecting the subspaces Vs = O^vl = and = 1, ^^l = 1 (which includes the 
first excited state) are 



2ns (2nL + 2L+1) 



2L + 1 



(26) 



{fis, Vs = 1; (fiL - 1) , J^L = 1| Tl^^ \ns, Vs = 0; n^, vl = 0) 
If we write each eigenstatc of the Hamiltonian as 



2hL i2hs + 1) 
2L + 1 



(27) 



(28) 



the matrix element of the T^^ operator between the ground state |0) = j'l'cOjO) and the first excited state |1) 
1*0,1,1) is 



2n,{2nL + 2L + l) q.o i,i 

C ' c 

Us ,nL Tls — I.'Hl 



2L + 1 



2hL {2ns + 1) 0.0 r.i 



(29) 



The formalism presented in this section provides the exact full solution of the problem. A simpler approach to 
study ground-state properties in the large N limit is provided by the mean-field analysis presented in the next section. 
This limit is a good benchmark to test more elaborate results. 

III. MEAN-FIELD ANALYSIS 

The geometrical interpretation of the Hamiltonian JQ) can be obtained by introducing a Hartree axial coherent 
state which allows to associate to it a geometrical shape in terms of a deformation variable (3. For a system with N 
bosons, this state is obtained by acting N times with a condensed boson F'l' on the boson vacuum |0) 



|7v,/3) = -=(rt)^|o) 



where the basic condensed boson operator is given by 

rt = 



1 



(30) 



(31) 



which depends on the f3 shape variable. The energy surface is defined as the expectation value of H in the intrinsic 
state 



EiN,p)^{N,P\H\N,P) 



1-x fl~l3^ 
X ^ „^ H — 



(32) 




FIG. 3: Energy surfaces per boson for the Hamiltonian Q for different values of the control parameter 



0.5. For X > Xc 
(broken phase) it is 



Minimizing the variational energy E(N ^ /3) with respect to (3 leads to a critical point at Xc 
(symmetric phase), the ground state is spherical and is obtained for f3 = whereas for x < Xc 
deformed since the minimum of the energy per boson eo — E{N, P)/N is obtained for /3 — y/l — 2x as can be seen in 
Fig. O At the critical point, it is worth noting that the energy surface is a flat surface near /3 = [2^l23l|. Within 
this mean-field (variational) approach, one thus gets the following ground-state energy per boson: 



ea{x > Xc) 
ea{x < Xc) 



l-x 



X 2 — 3x 



4 l-x 



(33) 
(34) 



One can also straightforwardly compute the expectation value of in the ground state which is found to vanish in 
the symmetric phase and equals 



(riL) = N 



1 - 2a; 
2(1 -x)' 



(35) 



in the broken one. 

However, other properties, as excitation energies or transition probabilities, that imply excited states require to go 
one step beyond this simple mean-field level. In the following, we shall use a combination of several methods already 
detailed for the simple case L = in Ref. which allow us to compute the corrections to these mean-field results 
as well as the gap or the transition rates which require the knowledge of excited states. 



IV. THE SYMMETRIC PHASE (1/2 < a; < 1) 



The starting point of our analysis is the elimination of the s boson by means of a Holstein-Primakoff boson expansion 
psf of s and L bosons (for a review in boson expansion techniques see Ref. [l^). Therefore, we introduce a set of 
bosons such that the mapping 



Lis 



= bib.. 

iVV25t(l-n,/iV)V2^(s%)t, 

= N-m. 



(36) 
(37) 
(38) 
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fulfils the commutation relations at each order in N in the Taylor expansion of the square root. 
With these notations, we have: 









pIPl 


= P^Pb, 




PlPs 




{N~2nb)+: nf :, 


pIPs 








= ^E(- 


-l)'^6t(l-n,/iV)W^(l-n,/7V)i/2 



(39) 
(40) 
(41) 



(42) 



where : A : denotes the normal-ordered form of the operator A. 

The Holstein-Primakoff mapping eliminates the s boson at the cost of introducing infinitely many boson terms. 
However, each term in the expansion has a definite order. As shown in the preceding section, for 1/2 < a: < 1, 
the number of L boson in the ground state goes to zero in the thermodynamical limit. Thus, to capture the finite 
N corrections, one performs a 1/A^ expansion of the Hamiltonian considering that ni,/N <^ 1. In this phase, the 
Hamiltonian Q reads: 



H 



1 - a; 
4 

1 - X 
4 

1-x 
4 



2(3x-l) / t 



1-x 



PlPb - 



-PlPb^^l^[Pl+Pb 



Pi 



+ Plm + ubPb 
+ Plnb + mPb 



0{l/N^). 



(43) 



Here, we have restricted this expansion to order (1/NY but the method we used can, in principle, be applied beyond 
this limit as shown in Ref. 2413 for the LMG model (L = 0). Our aim is to diagonalize H order by order. At leading 
order, one obviously recovers the mean-field ground-state energy per boson eo(A) = -f 0(1/A). At order (l/A)'^, 
the Hamiltonian is quadratic and it can thus be easily diagonalized_through a Bogoliubov transform giving rise to 
the boson Random Phase Approximation formalism presented in |27j and more recently exploited to describe the 
properties of the symmetric and broken-symmetry phases in the interacting boson model '21!, l28l| . Higher-order terms 
cannot be diagonalized by a Bogoliubov transformation, so that one has to resort a more sophisticated method. 



A. CUTs formalism 

The Continuous Unitary Transformations (CUTs) technique has been conjointly proposed by Wegner [29| and 
Glazek and Wilson [sollsH . For a pedagogical introduction, we refer the reader to Refs. 32, fs^l- Here we only sketch 
the main lines of this simple and powerful approach. 

The idea of the CUTs is to diagonalize the Hamiltonian in a continuous way starting from the original (bare) 
Hamiltonian H = H{1 = {)). A flowing Hamiltonian is thus defined by 

H{l)^UHl)H{0)U{l), (44) 

where I is a scaling parameter such that H{1 = oo) is diagonal, and U{1) is a unitary transformation, i.e. satisfying 
U{l)W{l) ~ U^l)U{l) — 1. Taking the derivative of Ea. (|H|l with respect to I yields the differential (flow) equation 

diH{l) = [,jil),H{l)], (45) 

where the generator of the unitary transformation ri[l) is 

r,{l) = diU\l)U{l) = -U\l)diU{l). (46) 

CUTs are also a powerful tool to compute the expectation value of any observable VI. As for the Hamiltonian, we 
define a flowing operator 



= U\l)VL{Q)U{l), 



(47) 
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which obeys 

din{i)^[v{i),n{i)], (48) 

with n ~ = 0). The expectation value of on an eigenstate of H is then given by: 

{tP\n\ij) = {i,\U{l = oo) = c3o) U\l = oo)|i/;), (49) 
where U\l — oo)|-!/)) is simply the eigenstate of the diagonal Hamiltonian H{1 = oo). 

The keypoint of this approach is an appropriate choice of the generator 77 which, in fact, depends on the problem 
under consideration. Here, the Hamiltonian H expressed in terms of h boson can be schematically written as: 

iJ(0) = i7o(0) + H+{{]) + H-{Q) + H+{{)) + ^2^(0), (50) 

where = (^1^2)^ 0, 1 or 2 subscripts indicate the number of created (+) or annihilated (— ) excitations. 

To perform the CUTs, we choose the so-called quasi-particle conserving generator first proposed by Mielke [s^ in 
the context of finite matrices and generalized to many-body problems by Knetter and Uhrig 35] which reads 

7,(0 = H+{1) - H^{1) + H+{1) - H^{1). (51) 

In the symmetric phase {Hf^ = 0) this choice coincides with the generator proposed by Stein [s^ . The flow equations 
are then simple quadratic functions of the Hamiltonians: 



diHoil) = 2[ [HUl),H-il)] + [H+{l),H-{l)]y (52) 
m+il) = [H+il),Ho{l)]+2[H+il),H^il)], (53) 
diH+{l) = [H+il),Hoil)]. (54) 

In the limit I — 00, the Hamiltonian conserves the number of b boson so that Hf2i^) = ^-iid H{oo) = Ho{oo). 
Following the method developed for the LMG model in '^^ convert these equations, which deal with operators, 

into equations involving coupling constants. This is achieved by expanding Hamiltonians Hq and powers of 



B. Flow equations for the Hamiltonian 



(see Sec. WV^ 



In the symmetric phase {Hf^ = 0), we have three elementary operators : nt, :,P^,Pi from which Hq and H2 are 
built. More precisely, the expansion of these Hamiltonians can be written as: 

^0{l) = 2^ ' (55) 

^ h^^^l Jl)Plpf : n? : 
HUD - E ' (56) 

a,/3,i5eIN 

Note that for L = 0, one has P/Pfc =: :, so that h^l\ ^ = ^k^a+2f3 o- ^^'^ then readily recovers expressions given 
in Ref. for the case of a scalar b boson. Using this expansion and Eas. H52l54|) . we can easily derive the flow 
equations for the couplings li'^^ ^{l) which are given in Appendix fXl up to order [l/N)"^. These flow equations can be 
solved exactly and, at order (l/A^)^, one finally gets: 

H{go) = /lo. 0,0(00) -I- /lo. 1,0(00)^6 -f /io,2.o(oo) : nl : +ho.o,i{oo)P^ Pb + /io,3.o(oo) : njj : -f /io,i,i(oo)Pjnf,P6 + 0{1/N^) 

(57) 

with 

hoi (0 

SeTN 

and 



4ao(oo) 
4ao(oo) 



1-a; 
2L+ 1 



1 - 3x 



-(2L + 5) + (6L + 13)x L + 2 



16S(x) 



4S(a;)i/2 



(2i + l)(l-a;)a; 

-{2L + - x)x^ 

{2L + 1) - (8L2 + 6L - 33) x + (32^2 - 2L - 149) a;^ - (24^2 _ 38L - 179) ; 

128S(x)5/2 

2L + 5 + (2L2 - 3i - 17) a; - (21,2 _ 5^ _ 20) ; 
16S(a;)2 



<l_o(oo) = S(a.)V2, 
4!i,o(oo) 



(1 — a:;)a; 



-l + (2L + 5)a; L + 2 



4S(a;) 2S(a;)i/2 



(1 - x)x^ 



L + 1 - (2L2 + 3L - 6) X + (12L2 + 15L - 23) x^ - (10L2 + 5L - 32) ; 

16S(x)5/2 

1 + (2L2 + 5L - l)a; - (2L2 + 3L - 4) a;2 
4S(^j2 



(0) _ (1 - x)x^ 



holfii'^) = -{l-x)x^ 



1 - 3a; + (12L + 29)x2 - 3(4L + 9)a;3 L + 1 - La; 



325(^)5/2 



— X 



4S(a;)2 



a;(l-.r)(:lr-l) 



8S(x) 

4!o,iM = -{'^-x)x^ 



l-3x + llx2-9a;3 1 + (L - 3)a; - (L - 4)a;2 

{2i1j ~\~ O) ~ r In 



1285(x)5/2 



8S(:r)2 



4°3,o(oo) 



{l-xfx^ 
8S(a;)5/2 ' 



(l-a;)2a;2 (l - 2a; + 9^2) 



64S(x)5/2 



where we have set ^(a;) = a; (2a; — 1). 

The ground-state energy per particle is thus given by 



eo(A^) = C,o(oo) + ^4!^,o(o«) + ^Co(o«) + ]^Co(o«) + 0{l/N% 
whereas the gap reads 

A(iV) = 4°l,o(oo) + ^<io(oo) + ^41o(oo) + 0{l/N% 
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Of course, these expressions coincide for L = with those given in Refs. MM- 

For L = 2, one recovers the results 
given in Ref. [s^l- The mean-field result ^6'6\ is also recovered in the thermodynamical limit. 

It is important to note that the Hamiltonian H{oo) — Hq{oo) is not diagonal in the eigenbasis of rib (except 
for L — 0) even though it always commutes with ni,. Consequently, for each number of excitations, H must be 
diagonalized. 

As can be observed in Eqs. H59(l - (|71|l . some divergences appears, for x = Xc, in the sub-leading corrections. We will 
see in Sec. IVII that the structure of this singular 1/A^ expansion at the critical point allows us to extract nontrivial 
scaling exponents whose determination is one of the main motivation of this work. 

C. Flow equations for 6|, 

We now proceed to derive the flow equation for the operator bj^{l) from which any other observable can be obtained. 
Analogously to the treatment of the Hamiltonian flow equations, the first step consists in transforming the flow 
equation (|48|l for = bj^{l) into a set of flow equations for couplings. Therefore, we expand bj^{l) in power of 1/A^. 

Generically, one may expect to generate any terms bj^"" P^^ : : Ph^b^ . Here, we shall restrict our discussion to 
order 1/iV for which one only has eight operators 

bl{l)=A^[l)bl + A^{^)b^ + B^{l)blnb^B^{l)n;b,^^ 

(74) 

where we have introduced fe^ = (— 1)^6_^ and where, as previously, each fmiction has a canonical 1/A^ expansion 
given by the number of bosonic operators it is associated with, namely 



A±(/) 
D±{1) 



0(l/iV2), 





(0 


N 






(0 


N 






(0 


N 



0(l/iV2), 
-0{\/N^). 



(75) 
(76) 
(77) 
(78) 



The initial condition is of course given by: = 0) = so that one only has a nonvanishing initial coupling which 

is A^^\q) — 1. The flow equations are then obtained for these couplings order by order using Ea. (|^ . The full set 
of equations is given in Appendix^ As for the couplings deflning the running Hamiltonians, these equations can be 
solved exactly and lead to 



4')(oo) 



2$(a;)i/4 
_(l-x) 
16a; 



2L + 3 2{L + 2) 



$(a;)V4 $(x)5/4 



(79) 
(80) 



Af{oo) 
A«(^) 



2 

16a; 



2L-F3 2{L + 2) 



$(a;)5/4 $(2:)3/4 



(81) 
(82) 



B(«)(oo) 



1 - X 



8a;$(x)7/4^ 
1-a; 
8x-I>(a;)5/4 ' 



(83) 
(84) 
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= 32.x^.^(x)V4 ^ (87) 

where we have set $(a;) = ^2^, = + F_) and F<j = - for each function F ^ A,B, C, D. 

The above expansion of 6jt(oo) allows us to compute for any operator il which can be expressed in terms 

of 6jj and for any eigenstates \^) and \^') of the Hamiltonian H{oo). In the following, we shall consider two different 
examples to show the power of this approach. 



D. Expectation value of the occupation number ul 

Let us first consider the case where ft = ul and where \'^) = \'i>') is the ground state of H. This quantity normalized 
by the number of bosons can be computed straightforwardly by means of the Hellmann-Feynman theorem which states 



^=a^[(l + 2/)Vo,o] 



(89) 



where y = x/{l — x). Since we have the expansion of /iq.o.o up to order (l/N)^, one can easily get (til) at this order. 
Here, instead, we compute it in terms of the flow equation for the operator bj^ obtained in the preceding section, using 
the fact that: 



(oo) = mioo) = ^6j,(oo)6^(oo), 



(90) 



where &Ji(oo) is given by Ea. l74() with final values Eas. (|79() - (|88|l . The ground state of the Hamiltonian (|57|) being 
defined as the zero b boson state |0), one has : 



N 



1 Y.AUoomb,bl\0)+Oil/N'), 



N 
1 



(2L + 1) aL"'(oo)^ + ^A^^^(oo)AL'^(oo) +0{1/N'^) 



iO)i 



N 

(2L + 1) 

N 



N 



3x - 1 
4S(x)i/2 



(2L + 1) x{l-xy' 



iV2 
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{2L + 3)x 2{L + 2) 
S(x)2 + S(a;)3/2 



Oil/N'), (91) 



where, as previously, ^{x) = x{2x — 1). It can be easily verified that this expression coincides with Ea. H89|) . 



E. Transition probability between the ground state and the first excited state 

As explained above, the real power of the CUTs method is that it allows to easily compute off-diagonal matrix 
elements of any operator between any eigenstates of the Hamiltonian provided one knows the expression of the 

I 1 2 

associated running operator. As an example, we focus here on the transition probability T — \{l\T[^g\0)\ between 
the ground state |0) and the first excited state |1). The operator T^^ was defined in Eq. It is important to note 

that here, the ground state has a zero angular momentum whereas the first excited state has an angular momentum 
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L. To determine the matrix element of Tl„ of interest, we shall proceed as for the occupation number and consider 
its 1 /N expansion in terms of the b boson: 



Tr 



(92) 



bl + bo 



2N 



blnb + mbo) +0{1/N^) 



To be consistent, given that we only have the expression of &J(oo) + 6o(oo) at order 1/N, we need to consider 
6j(oo)n(,(oo) + ni,{oo)bo{oo) at order (1/A^)''. Using the expression H74() of the operator bj, one easily gets 



(l|6t (oo) + 6o(oo)|0) = A+(c») + A_(oo), 



(93) 



and 



(l|6j(oo)n6(oo) + nb(oo)6o(oo)|0) = {2L + 3)A+{oo)A^{oof + A_(oo) [{2L + 2)A^{oof + A+{oo)'^] . (94) 

Truncating these expressions at order 1/N and respectively and using Ea. H75|l and Eas. H79|l - H8Q(l . one finally 

obtains: 



T 

N 



S(x)l/2 



X 



(2L + 1) - 4(2L + l)x + (lOL + 7)x^ ~L + (3L + 2)x 



2S(a;)3/2 



0{1/N^). 



(95) 



As for the expansion of the spectrum, some singularities appears, and we shall see that they also provide the scaling 
exponents at the critical point. 



V. THE BROKEN PHASE (0 < x < 1/2) 



As shown by the mean-field analysis, for x < 1/2, the order parameter {ni,)/N has a nonvanishing value. This 
implies that we have to consider a new vacuum for the Holstcin-Primakoff expansion. Therefore, we shift the bosonic 
modes by a term proportional to We thus define the c bosons by 



bl 



(96) 



where the A^'s are complex numbers which form a (2L + l)-dimcnsional vector. Of course, the symmetric phase 
results are recovered when setting = 0. Then, using definitions H3t)|) - (|38|l and assuming that ric/N ^ 1, we expand 
the Hamiltonian that now contains some a term proportional to ^/N which reads 



(c^-A + At.5) + [-2(l-nA) + Pi+PA]| + ^ [(c^-At) (P;, _ 1 + „;,) + (A • c) (p| - 1 + n^) 



(97) 

where = A;; and A^ = (-1)^A_^. There are several choices of the A^'s which allows one to get rid of these terms. 
Here, we have chosen to set A^ = X^d^fi with Aq = 2{i~x) ■ ^^^e that in the thermodynamical limit, we recover the 
mean-field value (|35|l 



N 



E 



1 - 2x 
2(1 -x)' 



(98) 



Further, we emphasize that this choice of the A^'s is the same as the one proposed in the mean- field analysis where 
we have broken the spherical symmetry by populating macroscopically fi ~ boson level only. With this choice, the 
Hamiltonian reads: 



H 



'Nx 



3x- 2 

4(1 -a;)" 



(1 — 3x)(l — 2.t) X 5, , t X, i s 3, ^. +2 2, 

^ g^/_^) ' + + ^(1 - 2x)clco - -{p! + Pc) + g(l - 2x){ci + Co) 



-N' , + ^n, + -fl - 2x)clcn - -(p1 + P r) + -(1 - 2x)(cr + Cn) +0(l/\/]V). 

(99) 

Contrary to the symmetric phase, we do restrict our discussion to this order because, as we shall see later, the existence 
of gapless modes at this level does not allow to go beyond with the CUTs. 
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A. The spectrum 

The Hamiltonian H99|l can be easily diagonalized via a Bogoliubov transform. Therefore, we introduce the d bosons 
defined by: 



= cosh(e^/2)4 + sinh(e^/2)J^, 
= sinh(e^/2)4 + cosh(e^/2)4. 



(100) 
(101) 



The angles 0^ are chosen so that H written in terms of the d's is diagonal. From Ea. (|99() . it is clear that modes with 
and = plays a different role and actually decouple. As can be easily seen, eliminating off-diagonal terms 
for n implies to set Qfj,:^o — oo and gives 2L gapless modes. Since such a transform is singular, one has to use 
another route. The contribution of terms with ii =/= in the Hamiltonian reads 



Introducing the position and momentum operator 



V2 



Pn 



V2 



one has 



X 

2 



(102) 



(103) 



(104) 



Since [P^ + (— 1)'^P_^, — (— 1)^X_^] = 0, i/^ is written in a diagonal form and its spectrum is indeed found to be 
gapless and continuous. The correction to the ground-state energy coming from this contribution is thus —Lx/2. 
Let us now consider the fi — part of the Hamiltonian which reads 



- 2x 



CoCo+ I --xl (c^ +Co). 



For this contribution, the Bogoliubov transform can be simply achieved and one gets: 



H, 



^Vl - 2a; - ^ + 2x^ + Vl - 2x djdo- 



(105) 



(106) 



The correction to the ground-state energy coming from this contribution is thus given by the do boson state. As a 
result, Co at this order reads: 



, , a; 2 — 3a; 1 

eoiN) = -- + - 

4 1 — a; A* 



-2-2(L-2)a;+(2L-l)a;2 1 



4(1 - a;) 



-Vl - 2a; 
2 



0{1/N^ 



(107) 



As in the symmetric phase, the leading corrections coincide with the mean-field result 1)3411 and L only appears in the 
sub-leading terms. 

Concerning the gap, the above analysis indicates the existence of 2L gapless modes and a gapped one with excitation 
energy 



A'(A^) = Vl-2a; + 0(l/A^). (108) 
As previously, one can simply obtain {ni)/N by replacing /io,o,o by eo{N) in Ea. (|89|l . and the result is 



{til) _ l-2x 1 
~/V~ ~ 2(1 - a;) ^ 2N 



L 



VI - 2a; X - 1 



0{1/N'^). 



(109) 



At this stage, one can understand the difhculty to go beyond this order in the presence of gapless modes. Indeed, 
computing the next-order corrections would imply to keep on making the 1 /N expansion around the (broken) vacuum, 
but such a procedure does not take into account the degeneracy due to the gapless modes. Note that for L = where no 
gapless modes emerge, we have been able to compute these corrections using CUTs 0. To conclude this subsection, 
we wish to underline that in the two-level BCS model where gapless modes also exist, Richardson has obtained the 
finite-size corrections in the broken phase beyond the Bogoliubov order using the 1 /N expansion of the exact solution 
[381] , whereas we computed them more recently using CUTs in the symmetric phase |39j| . 
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B. Transition probability between the ground state and the first excited state 

I 1 2 ' t 

As in the symmetric phase, we shall now compute the transition T — {1\Tl^ |0) where — Lo + L^s. However, 
the important difference is that, in the broken phase, one has 2L gapless modes which renders the definition of the 
first excited states more tricky. In the thermodynamical limit, the ground state thus becomes infinitely degenerate 
and one actually has to simply consider the expectation value of over the ground state. To avoid any confusion, 
we shall call this quantity T' instead of T. Using the expansion H93I) and the shift H96|) with the choice of the given 
previously, one gets 

r = |(o|Tijo>p, 

= nH\1{i-\1) + o{n), 

-^'t^^ + OW- (110) 
(1 — x)"^ 

Firstly, it is important to note that T' is proportional to iV^ in this phase whereas T scales as N in the symmetric 
phase. Secondly, in the broken phase, T' vanishes at the critical point whereas T diverges when approaching from the 
symmetric phase. This result clearly suggests an anomalous scaling behavior at the critical point that we shall now 
investigate in details. 



VI. THE CRITICAL POINT 



In this section, we shall analyze the behavior of the 1/A^ expansion of the quantities considered in this study: the 
ground state energy, the gap, the expectation value of tt-l in the ground state and the transition rate T between the 
ground state and the first excited state. The common point of all these expansions is that t hey become singular at 
the critical point. Following the arguments presented in a recent series of papers [l4L |24L 1371 l39j we shall now recall 
how this intriguing property allows one to extract the finite-size scaling exponents at this point. 

All quantities considered in this study display a singular behavior for x = Xc- This singular behavior can emerge 
in sub-leading corrections as for the ground-state energy but also in the leading term as illustrated by the transition 
rate in the symmetric phase [see Eq. (|95f) ]. Thus, schematically, the 1/A^ expansion of a physical quantity $ can be 
decomposed into a regular and a singular part as 

^Nix)^^'^H^) + <^7H^), (111) 

where, contrary to $5^^ and all its derivatives with respect to x do not diverge when x goes to Xc- Furthermore, 

at each order of the expansion, the divergence of is dominated by a single term. To be more concrete, let us 

consider the ground-state energy in the symmetric phase for which 



sing 
N 



1-x 1 (2L+ 1)(1 -Sx) 
~ ^ N 4 
1 (2L+l)S(x)V2 1 



'^Ti-) - ^ + -^— (112) 



^"^^ ^ N ^-^^ + ]^^o:o,o(«^) + ^/io:o.o(^) + 0(l/iV*). (113) 



In the vicinity of the critical point, these diverging terms have a leading contribution which is proportional to ^{x) ^ 
for /ig^Q_g(oo) and to S(x)~^/^ for /Iq^q q(oo). This has lead us to conjecture that near Xc the singular part behaves as: 



<^>7Hx) ^^:F^ NE{xfl^ , (114) 



where is a function that only depends on the scaling variable A^S(a;)^/^. We underline that for the LMG model 
{L = 0) we have checked this scaling hypothesis up to high order in the 1/N expansion. For the ground-state energy 
discussed above, one has f$ = 1/2 and n$ = 1. 

Once the form (|114|l is accepted, the scaling exponents are directly obtained. Indeed, since at finite N , physical 
quantities must not diverge, we conclude that, necessarily, Tij,{x) ^ 2;-244./3 gp ^-^^^ finally one has, <I>^"^(xc) ~ 
jY-(n4,+2^4,/3)_ TablelHwe have gathered the exponents obtained for the quantities discussed in this paper. 

We wish to emphasize that the scaling exponents related to the spectral quantities, i.e., eo, A and, using Eq. H89|l . 
(rii), can also be obtained in a different way. Indeed, as explained in Sec. IIIII the energy surface is the one of a 
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n$ 


-(na, + 25<i,/3) 


eo 


1/2 


1 


-4/3 


A 


1/2 





-1/3 


(til) 


-1/2 





1/3 


T 


-1/2 


-1 


4/3 



TABLE I: Scaling exponents for the ground-state energy per boson eo, the gap A, the number of L bosons in the ground state 
{ul) and the T transition probability. 

quartic oscillator (/3^-like potential) and this can be used as a starting point of a semi-classical description to show 
that the spectrum, at the critical point scales as N~^/^ . For technical details, we refer the reader to Ref. for the 
LMG model or |43| for the IBM with L = 2, and we also note that the "critical" scaling exponents do not depend 
on L. However, the present approach has a real advantage as compared to this latter method since one can compute 
the scaling exponents of any observables using the expression of &^(oo). We are further not restricted to expectation 
value but we can also investigate off-diagonal terms as illustrated with the transition rate T. 

VII. NUMERICAL RESULTS 

In this Section we check the validity of the analytical expressions obtained in the preceding Sections using CUTs. 
The observables studied are: the ground-state energy per boson eo, the gap A, the expectation value of the number 
of L bosons in the ground state (rti) and the transition probability between the ground state and the first excited 
state T. 

N=500, L=2 




0.5 1 0.5 1 

X 



FIG. 4: General features of the observables studied in this work as a function of the control parameter x obtained by numerical 
diagonalization. 

In Fig. 21 we present the general features of the selected observables as a function of the control parameter x for 
L = 2 and for N — 500. Note that in the broken phase (x < Xc), we have plotted A and T/N"^ instead of A' and 
T' /N"^, these two latter quantities being discussed in Sec. IVII Bl We can thus clearly appreciate the emergence of 
Goldstone modes in the broken phase. We emphasize that {nL)/N as well as the transition probability T/N"^ may 
be considered as order parameters since they vanish in the symmetric phase and acquire a finite value in the broken 
one. However, while {nL)/N is directly related to the physical ground state, T involves the first excited state which 
turns out to collapse into the ground state in the broken phase. This latter property makes it a more controversial 
candidate for an order parameter as recently discussed in Ref. p^l4ll| . 

In Fig. |31 we plot the difference between the numerical and the mean-field value of eg (dashed line) and (n^) /N (full 
line) as a function of the boson number, N , for three characteristic values of the control parameter: x = 0.75 in the 
symmetric phase, x — 0.5 at the critical point and x = 0.25 in the broken phase. As can be seen, the mean-field results 
become exact when increasing the number of bosons. It is interesting to note the change in sign in the deviations of 
the order parameter showing that the mean-field approach underestimates (overestimates) it in the symmetric phase 
(in the broken phase). To emphasize this effect, we plot in Fig. El the same quantities for = 20 bosons as a function 
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of the control parameter x. While deviations in the ground-state energy behave smoothly around the critical point, 
there is a well defined cusp in the deviations of the {ni) /N . 

0.04 
0.02 


-0.02 
-0.04 

-0.06q 20 40 60 80 20 40 60 80 20 40 60 80 
N N N 

FIG. 5: Differences between numerical (num) and mean-field (m.f.) results for the ground-state energy (per boson) eo and 
expectation value of the number of L = 2 bosons in the ground state (per boson) {nL)/N as a function of the boson number 
A'^ for three values of the control parameter: x = 0.25 (broken phase), x = 0.5 (critical point) and x — 0.75 (spherical phase). 
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FIG. 6: Differences between numerical (num) and mean-field (m.f.) results for the ground-state energy (per boson) eo and 
expectation value of the number of L = 2 bosons in the ground state (per boson) {til) /N as a function of the control parameter 
X at fixed = 20. 



Now that we have shown the main characteristics of the physical quantities of interest and the general agreement, 
in the thermodynamical limit, with the simple mean-field results presented in Sec. IIIII let us analyze in details the 
finite-size corrections in each phase independently. 



A. The symmetric phase 

In Section IV, we have obtained analytical expressions for the different corrections in the 1/A^ expansion of the 
selected observables. In order to check these results, we present several plots focusing in a first step, on the case L ~ 2 
and the dependence with x whereas, in a second step, we discuss the dependence with L. 

Let us first consider the ground-state energy per boson cq whose expansion in the symmetric phase is given in Eq. 
(|72|l . In Fig. the leading term in Eq. H72|l is compared with the numerical results for different N values confirming 

that the /iq^q q is indeed the true asymptotic value of cq. 

Next, we compare in Fig. ISjthe numerical and analytical subleading corrections to eg at each order. The numerical 

corrections of order p to cq are defined from the numerical value e™™ as e™™ — X]q=o ''■o^o o(°°)/-^" whereas 
the analytical correction is obviously given by /io^o_o(°°)- We present results tor p — 1,2,3 and N — 10, 100, 1000. As 
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can be clearly seen, for the largest value of = 1000, numerical and analytical results are almost indistinguishable 
even for values of x close to the critical point where /iq q o ^^"^ o o known to diverge. Note that the critical point 
Xc = 0.5 was explicitly excluded. Of course, the smaller N the larger the discrepancy since the numerical correction 
defined above still contains higher-order terms which play a role in this case. 

Along the same line, we analyze the corrections for the gap A defining the numerical correction of order p from the 

numerically calculated gap A""™ as N'' A""™ — Sa=-i ''■o"! o(°^) f-^" with il{oo) — 0. The analytical correction 
of order p is h'^l o(oo). 

Finally, we perform the same analysis for {ul) /N (see Fig. I10|l and T/N (see Fig. Ill|l by comparing the two first 
terms of their expansion with the numerical results. The numerical corrections are computed as for the gap. All these 
plots reflect that the x-dependence of the analytical expressions obtained with CUTs are in complete agreement with 
the exact numerical results for large values of N . 

To end up with these checks, we have investigated the L-dependence of the analytical results. We present in Fig. 
Elthe same observables (with the same notations) as those presented in Figs. ISIlll for fixed x = 0.6 and N = 1000 as 
a function of L. Once again, the agreement between the numerical results and the analytical expressions is excellent 
and confirms the validity of our analytical results. 
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FIG. 8: Comparison in log-normal scale between the numerical (symbols) and analytical results (lines) order by order for the 
ground-state energy per boson eo (see text for definitions). 
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FIG. 9: Comparison in log-normal scale between the numerical (symbols) and analytical results (lines) order by order for the 
gap A (see text for definitions). 
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FIG. 10: Comparison in log-normal scale between the numerical (symbols) and analytical results (lines) for the expectation 
value of the occupation number in the L level per boson in the ground state {nL)/N. n^^' stands for the 1/A'^ term and n"^^ 
stands for the term in Eq. ||MJ. 
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FIG. 11: Comparison in log-normal scale between the numerical (symbols) and analytical results (lines) for the transition 
probability per boson between the ground and the first excited state T/N. T^"^ stands for the A''-independent term and T'^-* 
stands for the 1/A'^ term in Eq. Il95t . 
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FIG. 12: Comparison in log-normal scale between numerical (symbols) and analytical results (lines) as a function of L. 
Notations are the same as in Figs. ISlTTl 



B. The broken phase 

As explained in Sec. the presence of Goldstone modes in the broken phase prevents from computing the 
corrections at high order. Therefore, we restrict our discussion here to the first nontrivial order. We present in Fig. 
1131 a direct comparison between numerical results (symbols) and the analytical ones given in Eqs. (|107I108I109I and 
I110|l (lines) as a function of the control parameter x for N — 1000 and L = 2. 



N=1000,L=2 




FIG. 13: Comparison between numerical (symbols) and analytical (lines) results. We only plot here the leading terms for each 
quantities. 

It is worth reminding that the gap associated with a one-phonon state in the symmetric phase turns into a Goldstone 
boson in the broken phase. The first excited state in the latter phase thus corresponds to a two-phonon state in the 
symmetric phase. However, this gapped mode ()108ll goes to zero at the critical point in the thermodynamic limit 
N = oo. As in the symmetric case, one can appreciate the agreement between analytics and numerics as already 
discussed, at this order, in Ref. [2ll |. 

We have also checked that the subleading terms of cq and (ni)/iV (beyond the mean-field results), which contains 
a nontrivial dependence with L, were fitting with numerics. In Fig. 1141 we show for L = 2, a comparison between 
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numerical and analytical results for N = 10, 100, 1000. As in the symmetric phase, the larger N the better the 
agreement. The dependence with L is tested in Fig. ^|at fixed x and N. 




FIG. 14: Comparison between numerical (symbols) and analytical (lines) results for L = 2. As in the symmetric phase, we 
have substracted from the numerical data the leading term given by the mean-field treatment. 
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FIG. 15: Comparison between numerical (symbols) and analytical (lines) results for the subleading corrections. As in the 
symmetric phase, we have substracted from the numerical data the leading term; c'q ^ and n^^' refers respectively to the term 
proportional to 1/A'^ in Eq. 11U7II and Eq. 110911 respectively. 



C. The critical point 



We now turn to the critical point study. To check the value of the finite-size scaling exponents derived in Sec. IVII 
we have performed diagonalizations for large system size (up to iV = 2^^ bosons). Let us recall that for L = 0, we 
have checked these values for larger system size in Ref. We show in Fig. [H our results for different values of 

L = 0, 1, 2, 3. Note that we plot the log2 of each quantity. 
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FIG. 16: Plot of the singular parts of eo, A, {ul) and T at the critical point Xc = 1/2 as a function of the boson number A'^ 
for different values of L. 

In this figure, only the singular part of the physical quantities of interest is plotted, the regular one being removed 
using the ad hoc expressions given in this work. As can be seen, the exponents are independent of L as expected from 
our calculations and match very well the predicted values. 

For L = 0, these exponents can also be obtained by noting that the LMG model can be seen as an Ising model in 
a transverse field with long-range interactions [42. . Then, the scaling variable N'E?/'^ is obtained from the upper 
critical dimension of the equivalent model with short-range interactions which is known to be dc = 3 in this case. For 
L 7^ 0, the two-level system studied in this paper cannot be simply mapped onto a short-range interaction model. 
Thus, it is rather a remarkable fact that the finite-size scaling exponents are independent of L. However, as explained 
in Sec. IVII this is due to the /3'*-like potential underlying the critical theory. 
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VIII. SUMMARY AND CONCLUSIONS 



In this paper, we have studied two-level boson models where the lower boson has a zero angular momentum (s 
boson), and the upper one, an angular momentum L. All these models are defined by the U(2L + 2) algebra, from 
which one can find chains of subalgebras going down to the 0(3) angular momentum algebra. When the Hamiltonian 
is written as a combination of Casimir operators of a chain of subalgebras, it is said that a dynamical symmetry occurs 
and the problem is analytically solvable. In this paper, we focused on the study of the quantum phase transition that 
appears when the boson system has a 0{2L + 1) symmetry, i.e. a transition from U{2L + 1) to 0{2L + 2) dynamical 
symmetries. This second-order transition is well described by a mean-field approach and the subtleties arise in the 
finite-size corrections. Here, we have explicitly computed these corrections for several physical quantities using firstly 
a 1/A'^ expansion naturally given by the Holstein-Primakoff representation of the angular momenta, and secondly 
the continuous unitary transformations to diagonalize the Hamiltonian. In the spherical (symmetric) phase, we have 
thus been able to capture corrections beyond the standard Random Phase Approximation and to show that the 1/iV 
expansion is singular at the critical point. The analysis of these singularities has allowed us to compute the finite-size 
scaling exponents which have been found to be independent of L. In the deformed (broken) phase, we have only 
computed the first corrections via a simple Bogoliubov transformations, in order to show the main difference with the 
spherical one. 

We have also presented a formalism based on boson seniority that provides a simple and efficient way of solving 
numerically the problem for a large number of bosons (a few thousands). Using this powerful algorithm, we have 
compared order by order analytical and numerical results and found an excellent agreement between both. We hope 
that the present work will help in understanding the approach to the macroscopic limit in such models, a problem 
that has recently drawn much attention [l^ lioj . 
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APPENDIX A: FLOW EQUATIONS FOR THE HAMILTONIAN IN THE SYMMETRIC PHASE 

In this appendix, we give at each order, the flow equations for the couphngs and the initial conditions obtained 
from the expansion of the Hamihonian H43|l . For clarity, we have not explicitly written the Z-dependence of all 

functions /i^*^^ ^. 

1. Order (l/TV)-^ 

At this order, one has one flow equation 

5/C.o = 0, (Al) 



with 



Co(0) = (A2) 



2. Order (l/N)" 

At this order, one has three flow equations 

dih^olo = -4(2L + l)4'?(l,o': (A3) 

dih^olo = (A4) 

'-'I'hflfi ^ ^''■o,i,o'''2,o,0' 1^3; 

with 

Co(0) = 0, (A6) 

<lo(0) - (A7) 

41o(0) = (A8) 



3. Order (l/A^)^ 



At this order, one has six flow equations 
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At this order, one has ten flow equations 
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with 
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APPENDIX B: FLOW EQUATIONS FOR bl IN THE SYMMETRIC PHASE 

In this appendix, we give at each order, the flow equations for the coupUngs involved in the 1/N expansion of bjj^{l) 
[see Eg . (1741) ] and the corresponding initial conditions. For clarity, we have not explicitly written the /-dependence 
of all functions. Further, it is convenient to introduce Fg — ^{F^ + F_) and F^ — ^{F^ — F_), for each function 
F = A,B,C,D. 



1. Order (l/iV)° 



At this order, one has two flow equations 



with 



4°nO)=4'"(0) = l/2. 



(Bl) 
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(B3) 



2. Order (1/iV)^ 

At this order, one has eight flow equations which decouples in two sets of four equations. 
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APPENDIX C: SOLVING THE FLOW EQUATIONS 

The flow equations given in the above appendices have to be solved order by order in At order (l/iV)~^, 

nothing has to be done, so let us turn to order (1/A^)°. The equations for h'j^l q(/) and h^l q(/) are easily solved by 
noticing that /iq°] o(0^ ~ '^^2^0 o(0^ ^ constant of the flow. One gets the hyperbolic solutions 



<l,o(0 

"■2,0,0 



(0 



tanh [2Aoo(/ + ;o)]' 
-sgn(e)Aoo 



(CI) 
(C2) 



2sinh [2Aoo(^ + /o)]' 

We have denoted Aqo the gap of the system at the thermodynamical limit, that is /iq^J q(oo) (see Eq. (j53J). The 
quantity Zq is such that the initial conditions are fulfilled, namely q(0) — Aqo/ tanh [2AooZo], and we also introduced 

the best way to solve the flow equations is in fact to 

(C3) 



h2'o_o{0)/[2h'^\ o(0)]- already explained in Ref. 



introduce a new "time scale" which is more adapted to the problem, and defined by 

t = sgn(e) exp [2Aoo(/ + lo)] , 
with initial conditions now given at 

to = sgn(e)exp(2AooZo)- 
After some algebra to can also be shown to be equal to 



1 



to = - i + vT 



Eqs. HC1|I and (|C2|) now read 
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h^°^ (t) - -A 

"2,0,0v''/ ^ 
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(C4) 
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(C6) 
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00 

This 



which are rational expressions in t. The renormalized values at Z — > cx) are now found by taking the limit t 
sgn(e)oo. Let us remark that the off-diagonal coupling h'^l Q{t) goes to zero and behaves like for t —i- t _ 
will be true for all off-diagonal couplings creating two excitations since the energy cost of such excitations, in the 
thermodynamic limit and for large t, is nothing but 2Aoo, so that the couplings must vanish as exp(— 2AooZ). 

The last flow equation (for the spectrum) at order (1/A^)° is solved by noticing that 2h'^l git) — (2L + l)hl^\ ^{t) is 
a constant of the flow, equal to its initial value, namely —{2L + l){3x — l)/2. 

The last task at order {1/N)° is to obtain the solution for the observable. For this, one simply has to insert Eq. (|C7|) 
in ljBl|l and IIB2p . then replace di with 2Aootdt, and solve the resulting equation. This yields 



1 / (t-l)(to + l) 

2V (i-Hl)(io-l)' 

1 l {t + l){to-l) 
2V {t-l)ito + l)' 



(C8) 
(C9) 
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The next orders are solved in the same fashion : one inserts the expressions known from the previous orders, replace 
di with 2Aootdt and solve the equations (which is most simply achieved thanks to a computer algebra program). We 
refer the interested reader to the details given in Ref. j2^, where some more technical details are given. 

Let us emphasize that the usefulness of the t variable comes from the fact that there exists only one basic energy 
scale in the problem, namely the gap. All energy scales are integer multiples of this gap. As previously mentioned, 
off-diagonal couplings associated to an energy scale 2Aoo decay as t~^, and in the general case, an off-diagonal coupling 
whose energy scale is nAoo decay as Such a time variable would be useless in a problem where many different 

energy-scales exist. 
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